Ferromagnetic zigzag chains and properties of the charge 
ordered perovskite manganites 
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Abstract 

The low-temperature properties of the so-called "charge ordered" state in 50% doped perovskite 
manganites are described from the viewpoint of the magnetic spin ordering. In these systems, the 
zigzag antiferromagnctic ordering, combined with the double-exchange physics, effectively divides 
the whole sample into the one-dimensional ferromagnetic zigzag chains and results in the anisotropy 
of electronic properties. The electronic structure of one such chain is described by an effective 3x3 
Hamiltonian in the basis of Mn{3deg) orbitals. We treat this problem analytically and consider 
the following properties: (i) the nearest-neighbor magnetic interactions; (ii) the distribution of the 
Mn(3deg) and Mn(4^j) states near the Fermi level, and their contribution to the optical conductivity 
and the resonant x-ray scattering near the Mn iC-absorption edge. We argue that the anisotropy 
of magnetic interactions in the double-exchange limit, combined with the isotropic superexchange 
interactions, readily explains both the local and the global stability of the zigzag antiferromagnctic 
state. The two-fold degeneracy of eg levels plays a very important role in the problem and explains 
the insulating behavior of the zigzag chain, as well as the appearance of the orbital ordering in the 
double-exchange model. Importantly, however, the charge ordering itself is expected to play only a 
minor role and is incompatible with the ferromagnetic coupling within the chain. We also discuss 
possible effects of the Jahn- Teller distortion and compare the tight-binding picture with results of 
band structure calculations in the local-spin-density approximation. 

PACS numbers: PACS numbers: 75.25.+Z, 71.45.Lr, 78.20.Bh, 75.30.Vn 
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I. INTRODUCTION 



It appears that many perovskite manganites, when doped by holes up to the level cor- 
responding to the formal valence of Mn +3.5, exhibit very specific properties, inherent to 
only this particular hole concentration and different from other doping regimes. The com- 
pounds are insulating and at certain temperature range (typically at low temperatures) form 
rather peculiar magnetic ordering consisting of ferromagnetic (FM) zigzag chains, coupled 
antiferromagnetically - see Fig. Q In the case of three-dimensional manganites, this is the 
well known CE-type antiferromagnetic (AFM) ordering, which was proposed long ago by 
Goodenough,>i and first observed by Wollan and Koehler in Lai_a;Caa,Mn03 near x=0.5.— 
The zigzag AFM ordering appeared to be a more generic and was observed recently in a 
number of systems, including the layered manganite Lai/2Sr3/2Mn04 (Refs. OU) and the 
double-layer manganite LaSr2Mn207 (Ref. Q). The magnetic ordering is accompanied by 
the orbital ordering,=4 which is also shown in Fig. ^ The insulating zigzag AFM state can 
be transformed into the metallic FM state by applying an external magnetic field.— Other 
examples of such compounds and detailed description of their properties can be found in 
the review articles. - 

The systems are usually referred to as the "charge ordered" manganites, implying that 
the homogeneous state of Mn^'^"*" ions in the lattice is unstable towards the charge dispro- 
portionation amongst two alternately aligned types of sites: 2Mn'^'^"''— s^Mn^^+Mn^^. The 
electronic configurations of the Mn"^"^ and Mn'^"'" ions are and t^gCg, respectively. Thus, 
the orbitals are half-filled in both cases, giving rise to the AFM superexchange inter- 
action between nearest neighbors. Due to the strong Hund's rule coubling, the Cg electron 
and three t2g electrons at the Mn^"*" site are aligned ferromagnetically. It is assumed that 
the Cg electron is mobile and can hop to the neighboring Mn^+ by forming the FM double- 
exchange (DE) coupling in a certain number of Mn^^-Mn^"*" bonds.— The occupied Cg orbitals 
are ordered in the direction of these bonds so to maximize the DE coupling. This is the 
basic physical idea underlying Goodenough's conjecture for the CE-type AFM ordering.— 
The picture had no solid justification and for many years remained to be on the level of a 
hypothesis. 

Only very recently it was realized by several authors that, at least at very low tempera- 
tures, the basic proporties of the " charge ordered" manganites can be understood exclusively 



3 



from the viewpoint of the DE physics.— liSiil Due to the strong Hund's rule couphng, the Mn 
ions form a high-spin state in manganites. Therefore, the intraatomic exchange sphtting is 
large. If it is much larger than the kinetic hoppings between nearest-neighbor Cg orbitals, the 
minority-spin Cg states can be projected out of the problem. In the case of AFM ordering, 
this will suppress all kinetic hoppings between Mn-sites with opposite directions of the spin 
magnetic moments. In the following we will refer to this picture as the DE limit. Thus, the 
CE phase in the DE limit will be effectively devided into the one- dimensional zigzag chains. 
This will create two geometrically different types of Mn-sites: the so-called bridge sites A 
and A' , and the corner sites B and B' - see Fig. ^ However, the difference is not necessarily 
related with the integer change of atomic charges when going from one type of Mn-sites 
to another. In the other words, the bridge sites are not necessarily occupied by the Mn^+ 
ions and the corner sites - by the Mn"^"*" ions. The two-fold degeneracy of eg levels plays a 
very important role in the problem. The isolated chain behaves as a band insulator. The 
only occupied Cg orbitals at the sites A and A' are Sx^— and 3?/^— r^, respectively, that 



is exactly the orbital ordering shown in Fig. ^ In our recent Letter (Ref. ilu ) we argued 
that the local stability of the zigzag AFM state can be explained by the anisotropy of the 
interatomic magnetic interactions in the DE limit, whereas the charge ordering itself may 
be irrelevant to the problem. 

In the present paper we will further ellaborate this DE picture. In Sec. |n] we will 
derive the tight-binding (TB) Hamiltonian for the zigzag chain, discuss basic features of the 
electronic structure, and obtain several useful expressions for matrix elements of the Green 
function and for the wavef unctions. In Sec. IIIII we will derive analytical expressions for 
magnetic interactions in the DE limit and address the problem of local and global stability 
of the zigzag AFM ordering. These two sections will also clarify the main results published 
m Ref. y. In Sec. IIVI we will turn to the optical properties of the zigzag chain and 
discuss possible implications for the resonant x-ray scattering near Mn i^'- absorption edge.— 
In Sec. El we will discuss validity of the TB picture by comparing it with results of first- 
principles band structure calculations in the local-spin-density approximation (LSDA). A 
brief summary will be given in Sec. IVIl 

We would like to emphasize from the very beginning that we treat the minimal model 
which can be proposed for perovskite manganites. The model is based on the DE physics and 
takes into account the two-fold degeneracy of the Cg levels. The combination of these two 
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ingredients represents a necessary basis for the analysis of manganite compounds. The main 
purspose of this paper is to argue that such a simple degenerate DE model provides a very 
consistent description for the low-temperature properties of the CE state in manganites. At 
present it is not clear whether the same DE model will be sufficient in the high-temperature 
regime. In Sec. IVH we will collect some remarks on this matter. 

We are planning this paper as the first part in the series of two publications dealing 
with the electronic structure and properties of the charge ordered manganites. The present 
work contains mainly the model analysis. Results of the first-principles band structure 
calculations and the structure optimization will be presented in Ref. The first reports 

nnn 

about this activity can be found in Refs. I13II14II151 In fact, the basic idea of the DE picture 
discussed in Ref. llOl was inspired by LSDA calculations of the zigzag AFM phase in layered 
manganites (Ref. |li 



II. ELECTRONIC STRUCTURE OF THE SINGLE ZIGZAG CHAIN 
A. Tight-binding Hamiltonian 

The basic translation properties of the single zigzag chain are given by the symmetry 
operation 5'i={C(7r) |Ro}, which combines the 180° rotation around the axis [1, 1,0], C{tt), 
with the translation Ro=[ao, oq, 0] connecting the sites A and A'. Therefore, it is convenient 
to choose the atomic bases at different sites in such a form that they can be transformed to 
each other by the operation S'„={C(n7r)|nRo}, where n is the integer number of translations 
separating two sites. For the Cg states this can be achieved by adopting the following basis 
orbit als: 

\1)a = -||3z2 - r')A + f\x'~ y')A = \Sx' - r')A 



(1) 



at the site A, 



at the site B, 



\l)B = \3z'-r')B 
\2)B = \x'-y')B 



\1)a' = -^13^2 - r')A' - f Ix^ - y')A' ^ \3y' - t')a' 
\2)a' = - r')A' + l\x'- y')A' = \x' - z')a> 



(2) 



(3) 
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at the site A', and 

\l)B' = \Sz'-r')B' 

\2)B' = -\x'-y')B' 
at the site B', which are periodically repeated to the whole chain. 

We assume that the kinetic hoppings between nearest-neighbor eg orbitals are of the dda 
type and can be parameterized in accordance with the Slater-Koster rules.— In the above 
basis, this yields the following nonvanishing matrix elements between sites A, B and A' of 
the chain 1 in Fig. ^ 

+11 - +11 - _i 

7-21 _ +21 _ 
''BA ~ ^BA' ~ 2 

In the notation t,^^,', n and n' stand for the site indices, L and L' - for the orbital (basis) 
indices. The absolute value of the dda two-center integral is used as the energy unit. The 
atomic orbitals |2) of the sites A and A' are inert and do not participate in the hoppings. 
In the local representation given by Eqs. (P)-©, any translation wRq in the real space will 
transform the matrix of kinetic hoppings to itself. Therefore, we can apply the generalized 
Bloch transformation C(n7r)e*"'^, — 7r<A;<7r, to the DE Hamiltonian in the real space 
if =||— This leads to the following 3x3 Hamiltonian in the reciprocal space, in the 
basis of generalized Bloch orbitals of \1)b{b') and \2)b{b') type: 



H{k) 



( Kl + e-^'^) -^(l-e-^'^)^ 

\{\^(^^) Ac 
y-^(l-e*'=) Ac } 



(6) 



The parameter Ac>0 was added in order to simulate the charge disproportionation between 
the sites Ai^A^ and B{B'\ The Hamiltonian has the following eigenvalues: 



e^=^(Ac±A/AQ + 8 — 4 cos A;), the bonding (— ) and antibonding (+) bands, and 
e°=Ac, a nonbonding band, which has nonvanishing weight only at the sites B and B' . 

All three bands are separated by energy gaps. Thus, if there is exactly one electron per 
two sites A and -B, the low-lying bonding band is fully occupied and the system behaves 
as an one-dimensional band insulator. Therefore, the insulating behavior of the so-called 
"charge ordered" manganites may not be related with the charge ordering. Rather, it may 
be a consequence of the peculiar zigzag AFM ordering, combined with the DE physics. 
Corresponding densities of states are shown in Fig. |21 (Details of calculations are given in 
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Appendix As expected for the one- dimensional system, the densities of states diverge 
at the band edges. Since only 3x^—r^ (3?/^— r^) state are occupied at the sites A{A'), the 
zigzag AFM ordering in the degenerate DE model automatically leads to the orbital ordering 
shown in Fig. 



B. Elements of the Green function 

In the reciprocal space, the Green function for the Hamiltonian (jHl) is given by 
G{k,e) = [e~H{k)+iri]~^, where 77 is a positive infinitesimal. It has the following matrix 
elements: 



GAAik,e) = \ - \ , (7) 



1 [ 4 £ 

l + cosfc r 1 1 1 3(1 — cos fc) 1 



2(4 -^fc) l4(^-^fe +^^) ^ki^-4 + iv)) 2(2-cosfc)e- Ac + iry' 

(8) 

3(1 — cos /c) r 1 1 "1 1 + cos/c 1 



2(4 -^fc) l4(^-^fc +^^) ^fc(^-4 + ^'?)J 2(2-cosA;)e- Ac + ir/' 

(9) 

1 + e-''' r 1 1 



GWk,s) = --— -—. —— , (10) 

2(4 -^fc) U-^fe £-4 + ^^ 

.12.. _^ V^{l-e~^')\ 1 1 



G^Lik, e) = ^^-^ <^ T > , (11) 

v/Sisinfc r 1 1 / 1 1 \ 1 



2(4 -^fe) l4(^-^fc +^^) ^fc(^-4 + ^^) \^k 4J^-^c + iv 

(12) 



The Green function elements in the real space are obtained by the Fourier transformation: 

G'Xni^) = ^ r dke^-'G'^i^:{k,e), 
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where r and r' correspond to either A or 5 in Eqs. (f7j)- (fT^ . This is the local representation 
for the Green function, meaning that the form of the basis orbitals L and L' varies from site 
to site as prescribed by Eqs. ©-(III). 
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C. Wave-functions 



We search the wave-functions of the Hamiltonian © in the form: 

' cose ^ 



^2(l+cos k) 



■ sin 6 cos ( 



(13) 



y -y/2{l-COsfc) 

By substituting this expression into the secular equation {H — e)|\I/)=0, the wavefunctions 
can be found as follows. 

For the bonding and antibonding states we obtain: cot 0-i-=— tan | and 



tan 0±=ei^ /^/2 — cos k. Corresponding wave-functions are given by 



1^^)^ = cose±|3a;2 -r^)^ 



(14) 



at the site A, and 



V 2 — cos k \ 2 2 



B 



(15) 



at the site B. 

For the nonbonding states £:°=Ac we have: tan0o=v^ tan | and 6q=tx/2. Then, the 
wave-functions will have nonvanishing elements only at the corner sites. They are given by 



l^o) 



1 



B 



V3 

\/2 — cos k I 2 



(1 + e'') tan ^13^2 - r^)^ + ^(1 - e"") tan^^ ^\x' - y^)B 



(16) 

The wave-functions at the sites N{B'^ are related with the ones at the site A{B^ by the 
transformation: 

|^)A'(B') = e''C'(7r)|^)^(B). (17) 
III. MAGNETIC INTERACTIONS AND STABILITY OF THE CE STATE 



We consider small nonuniform rotations of the spin magnetic moments on a discrete 
lattice, characterized by the angles If is the direction of the spin magnetic moment 

at the site % corresponding to an equilibrium, the new direction is given by ei=e°+<5ej, where 



(18) 



In the second order of {<5(^j}, the total energy change takes isotropic Heisenberg form (Ref. 
\v\ . see also Appendix |B|) : 



and the parameters { J^} can be expressed through the second derivatives of the total energy 
with respect to the angles {Sip^}.-^ In the AFM state we have e°=e^ (e^=— e°) if the sites 
i and j belong the same (different) sublattices. 



Following the work by de Gennes (Ref. |18|) we assume that the total energy consists 
of two contributions: i?t[{ej}]=i?£)[{ej}]+_E5[{ej}], where i?i5[{ej}] is the kinetic (double- 
exchange) energy of the Cg electrons, and Es[{ei}]=—^J^J2{ij) ' is the energy of AFM 
superexchange interactions between nearest neighbors (ij). In practice, the coupling J^< 
is associated with the localized spins, and may also account for some interactions amongst 
Cg electrons beyond the DE limit. In this model, the nearest-neighbor magnetic interaction 
is the sum of generally anisotropic DE coupling Jj^ (the second derivative of the kinetic 
energy of the Cg electrons) and the isotropic superexchange coupling: Jij=J^+J^ . 

Parameters {Jjij} depend on the magnetic state in which they are calculated and define 
conditions of the local stability of this magnetic state. The magnetic ordering is stabi- 
lized globally if it has the lowest energy amongst all possible locally stable magnetic states. 
In other words, this is the magnetic ground state of the system. Although the total en- 
ergy change near {e°} has the form of the Heisenberg model - Eq. (fTT?|l . it does not neces- 
sary mean that the same expression is applicable for the absolute value of the total energy 
in the point {e°}, i.e. generally it does not hold Et[{e^}]=—^ Joj{eQ ■ e°). One example 
illustrating such a violation is the DE model, where only interatomic DE interactions in the 
FM bonds contribute to the kinetic energy. The total energy in the DE model can still 
be expressed in terms of {Jy}, but the form of this expression will be different from the 
Heisenberg model. 

A. Intra-chain magnetic interactions 

In the DE limit, magnetic interactions between nearest neighbors located in the same 
chain, J^-^^i (s^e Fig. [T]for notations), can be calculated as:— 

JZb, = -^Im r de {Gl^^{e)t],\^^ + Gl^^{e)tl\^^} . (20) 
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Using matrix elements of the Green function, Eqs. (jlUp and and the kinetic hoppings, 
Eq. (jSj), we obtain: 



The interaction Ja^b^ always positive and responsible for the FM coupling within the 
chain. The numerical value for Ac=0 is J^_^^^ (0)^0.174. 

The charge disproportionation will reduce Ja^b^- Indeed, the former is defined as the 
difference of atomic populations at the sites A and S, An=nA—nB, which yields:— 



The dependence JaiBi versus An is shown in Fig. El As An increases, JaiBi decreases. 
Therefore, the charge ordering suppresses the FM DE coupling within the chain and plays 
a destructive role in the magnetic stability of the CE-type AFM ordering. 

B. Inter-chain magnetic interactions in the x-y plane 

In the DE limit, magnetic interactions between antiferromagnetically coupled atoms are 
given by the following expression: ^'■'i^'-' 



where Tr^ is the trace over the orbital indices. In Eq. (j22j) . the sites i and k belong to one 
zigzag chain, and the sites j and m belong to another zigzag chain, different from the first 
one. 

If i=Ai and j=i?2 for the geometry shown in Fig. ^ we will have the following six 
nonvanishing contributions to the exchange coupling (j22j: 




(21) 





(22) 



1' 
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The nonbonding x'^—z'^ {y^—z^) orbitals of A {A') sites are unoccupied and can be 
dropped in the analysis of magnetic interactions of the single zigzag chain. However, they 
contribute to the inter-chain interactions. The corresponding element of the Green func- 
tion is G^^^^^ (£:) = (£:— Ao+i?])^"'^, where Aq is the energy position of the atomic x^—z^ 
{y^—z"^) levels. The remaining elements of the Green function listed above are given 
by Eqs. ©-Q- Corresponding hopping matrix elements in the local coordinate frame 
are t^\^^=t^,^^=i, t]i^^=-t]^,^^=-t%^=^, and ti%,=-|. Then, noting that 
G B2B2{^)=G BiBii^) , we obtain the following expression: 

JZb2 = ^I^/" de{{G'l^^ie) + 3Gl\^^{e)) {Gl^.ie) + 3Gl^^{e)) + {Gl^.is) - SG 

(23) 

The numerical value of this integral for Ao=Ac=0 is J^^^^ (0)^0.106, which is significantly 
smaller than the value of the intra-chain integral JaiBiI^)- Such an anisotropy of magnetic 
interactions in the DE limit readily explain the local stability of the CE-type AFM ordering 
in the x-y plane. Indeed, by combining J^^Bi and JaiB2 with the isotropic superexchange 
interaction , we obtain the following condition of the local stability: 

JIb2 < \J'\ < JIb,, (24) 

i.e., we require the total coupling Jj^+J^ to be FM within the chain and AFM between the 
chains. Since -^a-^b2(0)<J^^5^(0), this condition can be easily satisfied. 

In the two-dimensional case, this is the only condition of the local stability of the zigzag 
AFM ordering. In the three-dimensional case, the inequality shall be combined with 
remaining two conditions of the local stability in the bonds A1—A3 and B1—B3 connecting 
neighboring x-y planes - see Fig. [T] 

C. Inter-chain magnetic interactions between neighboring x-y planes 

In the local basis of atomic orbitals (0)-®, the kinetic hoppings between sites Ai and 
A3, along the z direction, are tAiA3 = i tMA3=t%As = ^i and tf,A3 = |- Then, applying Eq. 
for i=Ai and j=As, and noting that GAiAi{£)=GAiA3{^), we obtain;^ 

JZa3 = ^I^/" i^MA^i^) + QGl^M Gl^^ie). (25) 
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The numerical value of this integral for Ao=Ac=0 is J^^^^ (0)^0.076. 

For the 5-sites, the only nonvanishing hopping along the z-direction is ^^^^^=1. Then, 
applying Eq. for i=Bi and j=B2, and noting that G BiBii^)=G B^Bsi^) , we obtain: 

JE,b, = ^1^/3 (^IbM' . (26) 

For Ao=Ac=0, Jb.b3 '^^^ ^e estimated as Jp-^ p^fO)— 0.117, which is different from the value 
reported in Ref. llj. The reason is the following. In Ref. (ly, we employed a different basis 
of atomic orbitals at the sites B and B'. In that basis, the Green function GbiBi{£) has 
off-diagonal elements with respect to the orbital indices, which have been neglected in the 
analysis of magnetic interactions in Ref. The new values of the exchange integrals take 
into account these terms. The correction appears to be small for the in-plane exchange 
interaction J^iB-zi^) modifies it only in the fourth digit after the comma. However, 
the change of the inter-plane exchange interaction J^^^Bsi^) appears to be more significant. 
With this new value of Jb^b-j.^!^)-! main statement of Ref. |lO| appears to be even stronger. 
If satisfies the condition 

the FM coupling is stabilized in the bond Ai-Bi, and the AFM cou- 
pling is stabilized in the bonds A1-B2, ^1-^3 and Bi-B^. Since 
'^BiB3(0)=max {7^*^52(0), Ji^^^3(0),J|'^53(0)}<Ji'^5^(0), appearance of the three- 
dimensional CE-type AFM ordering can be fully explained by the anisotropy of magnetic 

imit. We will call Eq. ()27j) the "hard" condition of the local 



interactions in the DE 
stability. Similar to Ref. 



lOl . we can also consider the CE-type AFM ordering in the regime: 



JZb, < \J'\ < JE.b,, (28) 



when the total coupling in the bond -B1--B3 is ferromagnetic. Even in this case, the CE-type 
AFM ordering may remain stable due to the joint effect of magnetic interactions in the 
bonds Ai-As, Ai-Bi and A1-B2. We will call Eq. the "soft" condition of the local 

stability. 
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D. Local and global stability of the zigzag antiferromagnetic ordering 



In Ref. |20| we argued that in the DE hmit the kinetic energy for a colhnear magnetic 
state M is related with the interatomic DE couphng parameter in the FM bonds, J^{M), as 
Ed(M)=-22^^J/^(M), where z^-^ is the number of FM bonds. The superexchange energy is 
given by Es{M)=—^J^{z-i-i—z-^i), where z-^i is the number of AFM bonds. This represents 
the connection between the parameters {Jij, J^} and the total energy of the DE model. 

Let us consider the three-dimensional case. For the hole concentration x=0.5 we can 
expect three possibilities: the FM ordering (M=F, Z||=6, and 2;|j^=0), the A- type AFM 
ordering (M=A, z^^=4, and z^i=2), and the CE-type AFM ordering (M=CE,Z|-f =2, and 
Z|j^=4). The chain-like C-type AFM ordering is unstable for x=0.5^ and can be excluded 
from the analysis. Corresponding total energies are given by: Et{F)=—12J^{F)—3J^, 
Ei(A)=-8Jf^(A)-J^, and Et{CE)=-4Jls^+J^ . 

Conditions of the local stability of the FM and A- type AFM states are \J'^\<J^{F) and 
J^(A)<|J'^|< J^(A), respectively, where J{j(A) is the nearest-neighbor exchange inter- 
action between antiferromagnetically coupled atoms in the DE limit. "Hard" and "soft" 
conditions of the local stability of the CE-type AFM ordering are given by Eqs. (j77j) and 
(OHj) . respectively. 

Numerical values of the exchange integrals for x=0.5 and Ac=Ao=0 are J^(F)~0.093, 
J|l(A)~0.115 and J||(A)^0.089 .^^i^°i^^ Corresponding phase diagram is shown in Fig. |^ 
where we plot the total energies for each magnetic state on the interval inside of which this 
magnetic state is locally stable. We also added the data for two canted spin states which 
may exist for x=0.5. These solutions of the DE model can be regarded as intermediate 
states between: 

(1) the FM and A- type AFM ordering (the canted spin state of the FA type); 

(2) the A- and G-type AFM ordering (the canted spin state of the AG type).— 

From Fig. |3] we can conclude that the CE-type AFM ordering may coexist only with 
the A-type AFM ordering and the canted spin ordering of the AG type. The FM ordering 
and the canted spin ordering of the FA type may develop only for smaller \J^\, when the 
CE-type AFM ordering is already unstable. If | J'^|>0.112 (the crossing point of CE and A 
total energy curves), the CE-type AFM state has the lowest energy amongst the states CE, 
A and AG, and is the magnetic ground state of the system. 
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In the two-dimensional case, the phase diagram is very simple: the total energy of the 
zigzag AFM state (z^|=2, z-]i=2), ii^t (zigzag) =—4J^^^^, shall be compared with the to- 
tal energy of the two-dimensional FM state {z^=A, z^i=Q), -Ej(ferro)=— 8J^(ferro)— 2J'^, 
where J-|^|(ferro)~0.115.^*^ Conditions of the local stability are: | J'^|< J^(ferro) for the two- 
dimensional FM ordering, and Eq. (j24|) for the zigzag AFM ordering. Two states coexist in 
the interval J^_^Q^<\J^\<JY^{ieiio). If | J'^|>0.112, the zigzag AFM state has lower energy 
than the FM state and vice versa. 

E. Simple model for the Jahn- Teller distortion around the bridge sites 

We assume that the Jahn- Teller distortion (JTD) results in the shift of atomic 
Ip'—z^ I'ix^—r'^ [x^—z^ I'iy^—r'^) levels of the A (A') sites by ±|Ajt. This is equivalent to 
the following choice of the parameters of our original Hamiltonian: Ao=Ajt and Ac=|Ajt, 
i.e. in addition to the local splitting at the A (A') sites, we shift all the states upwards by 
|Ajt. Obviously, the constant shift does not affect the magnetic interactions. 

Thus, the JTD in such a picture will be accompanied by a charge disproportionation, 
and suppress the FM DE coupling within the chain (see Fig. E)). This is, in principle, a 
negative consequence which lowers the upper boundary of the local stability of the zigzag 
AFM ordering. Nevertheless, for relatively small distortions, the change of J^^Bx 
expected to be small. The effect of Ajt on the interatomic interactions in the AFM bonds 
A1-B2, A1-A3 and B1-B3 is more dramatic. The results are summarized in Fig. El 

In order to evaluate the magnitude of the oxygen displacement which results in the split- 
ting Ajt, we consider the following procedure:— (i) we start with a general TB Hamiltonian 
on an undistorted cubic lattice, in the basis of Mn{3deg) and 0{2p) orbitals interacting via 
nearest-neighbor pdcr hoppings, the distance-dependence of which is simulated by the Harri- 
son low: (pda)(xd^'^^'^] (ii) we switch on the JTD of the MnOg octahedra around the bridge 
sites, the Q3 mode,-^ (it is assumed that the total volume of the MnOe octahedron remains 
to be constant); (iii) we project out the 0{2p) states. This results in the on-site splitting of Cg 
orbitals, Ajt— 285 (in units of the effective dda two-center integral for the undistorted cubic 
lattice), where the oxygen displacement is described by the parameter 6={dL—ds)/{dL+ds), 
with di (ds) being the long (short) Mn-0 bond length in the distorted MnOg octahedron. 

Then, the behavior of magnetic interactions shown in Fig. El suggests that even a modest 
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distortion 5=0.02 (Ajt— 0.56, see also Ref. 1271 ) may reduce the FM contribution to the 
AFM bonds A1-B2 and Bi-B^ by up to 20%, whereas the couphng in the FM bond Ai-Bi 
remains practically unchanged. Therefore, the JTD significantly lowers the low boundary 
of the local stability for the zigzag AFM ordering. As the result, the range of parameters 
\J^\ which stabilize the zigzag AFM ordering significantly widens. 



IV. OPTICAL PROPERTIES OF THE ZIGZAG CHAIN 

In order to evaluate the optical conductivity tensor a{uj) for the single zigzag chain we 
take the following steps. First, we include the Mn(4p) states in our model and evaluate 
their admixture into the eigenfunctions of Hamiltonian (jH)) using the simplest perturbation 
theory. We assume that the main perturbation is caused by the kinetic hoppings between 
nearest-neighbor Mn(4p) and Mn{3deg) orbitals.— The hoppings are of the pda type and 
can be parameterized in accordance with Slater-Koster rules.-^ It is further assumed that 
the energy splitting between the atomic 4p and Sdcg levels, A^p^e^, is much larger than the 
corresponding bandwidths, so that the band dispersion can be neglected in the perturbation 
theory expansion (this is certainly a very crude approximation, especially for the Ap band, the 
validity of which will be discussed in Sec. I^)- Thus, the mixture of the Mn(4p) and Mn{3deg) 
states will be controlled by the parameter a={pd(7)4p_ej A4p_eg, where {pda)^p_eg is the 
Ap-Sdcg two-center integral.'^ Then, by knowing distribution of the Mn(4p) and Mn{3deg) 
states we can find matrix elements of the interband optical transitions e^^e" and e^— 
and express them through the intra-atomic dipole Ap-Sdcg matrix elements. Finally, the 
optical conductivity can be calculated using the Kubo formula. Due to the one-dimensional 
character of the problem, the solution along this line can be carried out analytically. For 
simplicity we will consider only the case of homogeneous charge distribution (Ac=0). 

The simple TB picture considered in this section is of course very oversimplified and 
cannot provide a good quantitative description. However, it appears to be very useful in 
the analysis of symmetry properties and the structure of the conductivity tensor in terms of 
partial contributions from different sites of the system. 
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A. Distribution of Mn(4p) states 



First, we consider the admixture of the 4p-orbitals \x) and \y) into the bonding and 
antibonding states of the Hamiltonian (jU)). Taking in to account the exphcit form of the 
wavefunctions at the sites A and A' given by Eqs. (jl4j) and (jl7|) . and considering the nearest- 
neighbor pda hoppings from the sites A and A' to the site B of the chain 1 in Fig. ^ we 
obtain the following corrections to the wavefunctions at the site B, in the first order of a: 

|5^±)b = -acose^ {\x)b + e%)B} ■ (29) 

Starting with the sites B and B' and considering the hoppings onto the site A, we obtain: 

, ^ ^ , sin 9± sink , , , . 

V ^ — COS k 

Similar analysis applied to the nonbonding state |\E'o) yields: 

\5^o)a = V3ta^^^^L=\x)A. (31) 
V ^ — cos k 

Since |\E'o) has nonvanishing elements only at the sites of B and B' type and the kinetic 
3deg-4p hoppings are restricted by the nearest neighbors, the correction |5\E'o) may exist only 
at the (neighboring) sites of A type. 

Corresponding corrections at the sites A' and B' can be found using Eq. p7|l . 



B. Matrix elements of the optical transitions 

The nonvanishing intra-atomic dipole matrix elements obey the following 
rules>22i {3x'^—r'^\x\x) = {3y'^—r'^\y\y)=(3, {3z'^—r'^\x\x) = {3z'^—r'^\y\y)=—^P, and 
{x'^—y'^\x\x)=—{x'^—y'^\y\y) = ^l3. Then, using results of Sees. Ill (^1 and IIV Al for 
Ac=0, we obtain the following matrix elements: 

{^^\xmo)A = {^-\y\mrA' = (32) 

V ^ V 2 — cos k 
for the bonding- nonbonding transitions at the sites A and A'; 

[3 1 



8 V2 — cos k 
(33) 
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for the bonding-nonbonding transitions at the sites B and B 

, \*., = -in,r-i 

\/2 — COS k 



{6^-\x\^+)a+{'^-\x\6^+)a = + = -^(^f^^^^ (34) 



for the bonding- antibonding transitions at the sites A and A'; and 

(35) 

for the bonding- antibonding transitions at the sites B and B' . 
C. Optical conductivity 

The inter-band optical conductivity tensor can be obtained using the Kubo formula: 

a^yHoc-/ dA;(^i|7|^f)(^f|7'|^i)5(cu-4 + £^), (36) 
^ Jo 

where 7(7')= x, y; and e\ stand for the initial state i="— "; and stand for the 
final state f="0" or "+". We consider only site-diagonal contributions in Eq. ()36j) . The 
integral can be evaluated along the same line as for the local densities of states in Appendix 
1X1 Then, we obtain the following nonvanishing contributions to the conductivity tensor: 



^4.M = <M (37) 



at the sites A and A' where 



a^^^u;) oc \—^=^^]L= (38) 
c^V(3-^')(^'-l) 

for the bonding-nonbonding transitions (we drop the common multiplier a^/3^), and 



atiu;) oc ^V(12-.;2)(^2_4) (39) 
for the bonding- antibonding transitions; 

af,(^) = -agiuj) = -a^,{uj) = aj,(^) (40) 

and 

a^'(c.) = a^(u;), (41) 
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at the sites B and B', where 




3 1 



(42) 



4a;2^(3-cu2)(cu2-l) 



for the bonding-nonbonding transitions, and 




1 



(43) 



for the bonding-antibonding transitions. 

The results are shown in Fig. IHl The conductivity spectrum consists of two parts orig- 
inating from the bonding-nonbonding and bonding-antibonding transitions, and spreading 
within the intervals l<a;<v^ and 2<uj<2\^, respectively. Taking into account that a re- 
alistic estimate for the effective dda two-center integral is about 0.7 eV (for the undistorted 
lattice, without buckling of Mn-O-Mn bonds),— they roughly correspond to the energy 
ranges 0.7<cj<1.2 eV and 1.4<cij<2.4 eV. We note that the components of the conductiv- 
ity tensor associated with the bridge {A,A') and the corner {B,B') sites are characterized 
by different optical axes. The tensors ^^{uj) and a^' (u) are diagonal in the basis of cubic 
axes x=[l,0,0] and y=[0, 1,0] shown in Fig. ^ In this basis, each of the tensors has only 
one nonvanishing diagonal element, correspondingly (J^x{u) and ayy{uj), connected by the 
symmetry transformation ()37|1 . Therefore, the total conductivity a^{iu)+a^' (lu) behaves as 
an isotropic object in the plane x-y. Conversely, the tensors {uj)=a^' (uj) are diagonal in 
the basis of orthorhombic axes a=[l, 1, 0] and 6=[1, 1, 0], which are rotated with respect to 
cubic ones by 45°. In the orthorhombic frame, they have only one nonvanishing diagonal 
element (T^(c<j)=cr^'(co')=2(T^(co') in the direction of propagation of the zigzag chain (the 
6-direction). This gives rise to an anisotropy of the optical conductivity, which has been 
observed experimentally in Lai/2Sr3/2Mn04 by using a birefringence technique.-^ 

The energy range corresponding to the bonding-nonbonding transitions roughly agrees 
with the position of the mid-infrared peak observed in Lai/2Sr3/2Mn04 around 1.0-1.3 

15^30 jjjg]-^ intensity of experimental peak may be partly attributed to the (quasi-) one- 
dimensional character of the problem and divergence of the optical conductivity at the band 
edges. Another possibility is the lack of the inversion symmetry at the corner sites B and 
B'. In this case, the atomic Sdcg and 4p orbitals may belong to the same representation of 
the local symmetry group. This is an additional channel of mixing of the Sdcg and 4p states 
at the corner atoms, which contributes to probabilities of the optical transitions. 
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The JTD around the bridge sites can affect the conductivity spectrum indirectly, through 
the induced charge disproportionation, in a way similar to the intra-chain DE interaction 
considered in Sec. 1111 El Since the charge disproportionation Ac>0 increases the energy 
gap between the bonding and nonbonding bands and decrease the width of the former, we 
expect the upward shift of the optical absorption center and narrowing of the spectrum in 
the region of the bonding-nonbonding transitions. 

If the charge disproportionation does not take place (Ac=0), the optical conductivity is 
expected to be the smooth function in the region of bonding-antibonding transitions and 
reveals the same symmetry properties as for the bonding-nonbonding ones. 

D. Resonant x-ray scattering near the Mn K-absorption edge 

The resonant x-ray scattering near the Mn /^-absorption edge is considered as a powerful 
experimental tool, which allows to probe the charge and orbital distribution in the perovskite 
manganites.-^^ The interpretation of experimental data is, however, hampered by the fact 
that the direct dipole transitions are allowed to the 4p states but not directly to the 3d states 
in the case of i^'-absorption. Therefore, the understanding of how the 4p states interact with 
polarized Sdcg orbitals represents a special interest .^^2*^ Here we would like to show how the 
geometry of the zigzag chain imposes some constraints on the distribution of the Mn(4p) 
states, which should be reflected in the lowest energy part of the fT-absorption spectrum of 
manganites with the zigzag AFM ordering. 

In accordance with the TB picture (Fig. |21), the first unoccupied band which is involved in 
the transitions is the nonbonding band composed of the Sdcg orbitals of the corner atoms B 
and B'. If the hoppings between Sdcg and 4p orbitals are restricted by nearest neighbors, the 
4p states will have nonvanishing weight only at the bridge sites A and A' in the nonbonding 
part of the spectrum. This interaction, given by Eqs. (jHlj) and ()17j) . results in the ordering 
of \x)a and \y)A' orbitals. Since the Cg orbitals of the bridge atoms do not contribute to 
the nonbonding band, there is nothing to re-expand over the 4p orbitals of the neighboring 
(corner) atoms, and the corresponding weight of the 4p states at the corner atoms will be 
zero. 

Thus, in accordance with our TB picture, the lowest energy part of the i^'-absorption 
spectrum probes the \x)a and \y)A' states of the bridge atoms, which carry the information 
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about the distribution of the Cg states at the comer atoms. 

V. COMPARISON WITH LSDA BAND STRUCTURE CALCULATIONS 

The purpose of this section is to discuss vahdity of the TB model for the zigzag AFM 
state by comparing it with results of first-principles band structure calculations in the local- 
spin-density approximation (LSDA). LSDA and its extension in the form of the generalized 
gradient approximation (GGA) work reasonably well in mang 

anites.^2.2^«2L38 There are 

several unresolved problems, such as the relative stability of the Jahn- Teller distorted in- 
sulating A- type AFM state and the metallic FM state in LaMnOs: the standard band 
calculations account well for various magnetic properties of LaMnOa with the experimental 
crystal structure,— however the full optimization of the crystal structure substantially 
reduces the JTD and makes the A-type AFM state energetically less favorable than the FM 
state.-^ To say the truth, similar problem was encountered also in Lai/2Sr3/2Mn04, where 
after the structure optimization in GGA the FM state had lower energy than the zigzag 
AFM state.— However, it is also true that at present the problem has no feasible solution. 
For example, the LDA+?7 (or GGA+f/) approach only worsen the picture of magnetic in- 
teractions in LaMnOs.— This may be related with the fact that the relative position of the 
oxygen 2p and Mn(3(i) states (the so-called charge-transfer energy) is incorrect in LDA+f/, 
that leads to a large error in the calculated parameters of magnetic interactions.^ Thus, we 
believe that the LSDA provides the most reliable reference point for the analysis of magnetic 
properties of manganites. 

The quasi-one-dimensional TB model is a substantial simplification for realistic material, 
in several respect, (i) The DE limit is not realized in perovskite manganites. So, in the 
zigzag AFM state the chains are not fully isolated from each other, and the system is 
not strictly one-dimensional, (ii) The distribution of Mn{3deg) states is affected by other 
factors, which are not included explicitely in the TB model, particularly by the oxygen 2p 
and Mn{3dt2g) states.''^ (iii) The Mn(4p) states generally form a very broad band. In this 
respect, our attempt to simulate the optical properties by starting with atomic 4p-levels 
and considering their interations with the Cg bands in a perturbative manner was a crude 
approximation. It allows to catch the basic symmetry properties of the optical conductivity 
tensor, but certainly not all the details and the absolute values of the conductivity itself. 
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We use the ASA-LMTO method in the nearly orthogonal representation.— The Wigner- 
Seitz sphere radii for nonequivalent sites were chosen from the charge neutrality condition 
inside the spheres. The calculations have been performed for the phenomenological virtual- 
crystal alloy Yi/2Sr3/2Mn04 using the lattice parameters of the real two-dimensional per- 
ovskite Lai/2Sr3/2Mn04.^'^ In addition, we also investigated the effect of the JTD around 
the bridge sites. The latter is characterized by the parameter 6 introduced in Sec. IIIIEI 

The partial densities of Mn[3deg) and Mn(4p) states are shown in Fig. [7| They reveal 
many similarities to the TB picture. There are aslo some differences. We note the following. 
First, there is the clear sphtting of the majority-spin Sdcg states into the bonding band 
located around —0.5 eV (the case of 6=0), the nonbonding bands around 0.7 eV, and the 
antibonding band around 1.7 eV. Two nonbonding bands are composed of the Cg states of 
the corner {B and B') atoms and the y'^—z^ {x'^—z^) states of the bridge A [A') atoms. Even 
without the JTD, the system is insulating and the band gap is 0.3 eV. This value is however 
smaller than the estimate 0.7 eV expected from the tight-bonding analysis (the value of 
effective two-center dda integral). The difference is caused by the additional broadening 
of the bonding and nonbonding bands in the LSDA. The JTD increases the band-gap (for 
example, for 5=0.03 the band-gap is 0.7 eV). Second, the partial densities of states in the 
bonding part of the spectrum clearly shows the characteristic peaks at the band edges, which 
is the signature of (quasi-) one-dimensional behavior. The bandwidth is however larger 
than (-\/3— 1) x0.7~0.5 eV, expected from the TB analysis. The density of nonbonding 
and antibonding state is modified more significantly. Particularly, the nonbonding and 
antibonding parts of the spectrum overlap with the minority-spin t2g band (not shown in 
Fig. [7j). The minority-spin t2g states may interact with the majority-spin Cg states of the 
neighboring zigzag chains either directly, via next nearest-neighbor hoppings (note, that the 
nearest-neighbor hoppings between the t2g and eg orbitals are forbidden in the cubic lattice 
- Ref. Q), or indirectly, via the oxygen 2p band. As the result, the "nonbonding" states 
associated with the corner atoms B and B' are significantly broadened (to be compared 
with the 5-peak expected from the TB analysis - Fig. Conversely, the y'^—z^ and x'^—z'^ 
states of the bridge sites A and A' are only weakly bonded, resulting in the high peak of the 
density of states around 0.5 eV. 

Due to the hybridization, the distribution of the majority-spin Mn(4p) states around 
the Fermi level repeat the characteristic bonding-nonbonding-antibonding splitting of the 
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Mn(3(iep) states. The nonbonding part of the spectrum is composed mainly of the x {y) 
orbitals of the bridge atoms A {A') . The contribution of the Ap orbitals of the corner atoms 
into this region is significantly smaller. 

The corresponding optical conductivity is shown in Fig. |H1 Some details of calculations 



can be found in Refs. 
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39l . Here we only would like to mention that in order to speed 
up rather heavy calculations (note that taking into consideration the zigzag AFM ordering, 
the unit cell of Yi/2Sr3/2Mn04 consists of 56 atoms), the Brillouin zone integration in the 
calculations of optical conductivity was replaced by a summation with the phenomenological 
Lorentzian broadening of 0.136 eV. This causes some artificial broadening of the spectrum. 
Nevertheless, we believe that all important features are preserved in it, and were not washed 
out by the broadening. As the test, one can consider the off-diagonal element (Jxy{u) in the 
low-energy part of the spectrum. It clearly shows two peak around 0.4 eV and 0.8 eV (the 
case of 6=0), corresponding to two peaks of the (quasi-) one-dimensional density of states 
in the bonding region (Fig. Ej), in a close analogy with the TB picture (Fig. IH)). 

Surprisingly however that the shape of diagonal conductivity in the xy-plane is very 
different, contrary to our expectations based on the TB model. In the low-energy part of the 
spectrum it has only one large peak at 0.7 eV. The difference cannot be simply attributed 
to the broadening effects, because otherwise similar broadening would be expected also 
for the off-diagonal component of the conductivity tensor. Unfortunately, on the level of 
first-principles band structure calculations, it is rather difficult to decompose the optical 
conductivity unambiguously into partial contributions and to elucidate the origin of the 
spectral shape. Most probably, the difference between the LSDA calculations and the TB 
analysis is caused at once by several factors outlined in the beginning of this section. 

In the TB model, two contributions to the conductivity spectrum in the region of bonding- 
nonbonding transitions, cr^^^uj) and (j£,(c<j), are comparable (see Fig. The off-diagonal el- 
ement of the conductivity tensor, which is responsible for the anisotropy of optical properties, 
is related with the diagonal one by Eq. (jinjl . Therefore, (7^y{w) is expected to be of the same 
order of magnitude as the total diagonal conductivity a^^{u!)+a^^{u!). However, in a more 
realistic situation, when the system is not strictly one-dimensional, the optical anisotropy 
is significantly reduced. For example, in the LSDA we have only \axy{u!)\r^OAaxx{^^) in the 
region of bonding-nonbonding transitions. 

The JTD around the bridge sites results in the upward shift of the low-energy peak and 
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the simultaneous decrease of its intensity. Two experimental groups reported somewhat 
different results for the peak position at the low temperature: 1.0 eV (Ref. llSl ) and 1.3 
eV (Ref. I30!). Both are larger than 0.7 eV obtained in the LSDA calculations without the 
JTD. The JTD 5=0.03 yields the new peak position 1.0 eV, in agreement with one of the 
experimental reports. However, the obtained intensity 700 fi^^cm^^ at the peak maximum 
appears to be smaller than the experimental one, which is about 1100 i7^^cm~^ .^^i^° Another 
effect of the JTD is the redistribution of the in-plane conductivity around 2.0 eV, in the 
region of bonding-antibonding transition. 

As expected for the layered perovskite compound, the in-plane (axx) and out-of-plane 
(czz) components of the optical conductivity reveal a strong anisotropy, in a good agreement 
with the experimental finding for Lai/2Sr3/2Mn04.— 

Finally, we calculate nearest-neighbor magnetic interactions in Yi /?Sr3/9Mn04 as the 
function of the Jahn- Teller distortion S, using the same approach as in Ref. |35|. For small 
6 we obtain J^^Bi(^)— 43+3135 meV and JaiB2('^)— 12—6805 meV. Therefore, even without 
the JTD, the magnetic interactions in the zigzag AFM state show a very strong anisotropy 
<^Ai_Bi(0)> Jyii_B2(0)) in a good agreement with the TB analysis. In the LSDA, a small JTD 
5>0.02 becomes indispensable to stabilize the AFM coupling JaiB2{^)<^ between neighbor- 
ing zigzag chains. 



VI. SUMMARY 



Properties of the so-called "charge-ordered" manganites are typically regarded as one 
of the most striking examples of the strong coupling amongst the spin, charge and orbital 
degrees of freedom. The coupling, which is reflected in the number of transport, structural 
and magnetic phenomena, is also believed to be responsible for the peculiar form of the 
ordered state realized at low temperatures. 

In the present work we tried to elucidate the cause of such an unusual behavior and to 
understand it from the viewpoint of the magnetic spin ordering. We found that the zigzag 
form of the antiferromagnetic spin ordering, combined with the DE physics, already provides 
a very consistent description for various low-temperature properties of the "charge-ordered" 
manganites. In the DE limit, the one-dimensional FM zigzag chains become effectively sep- 
arated from each other. The one-dimensional character of the chains leads to the anisotropy 
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of electronic properties. The anisotropy of nearest-neighbor magnetic interactions in the DE 
hnait, combined with the isotropic superexchange interaction , readily explains the local 
stability of the zigzag AFM state. Formally, the basic DE picture is sufficient to understand 
the local stability. The JTD significantly reduces the FM DE contribution to the AFM 
bonds, and in this way widens the region of local stability of the zigzag AFM ordering. It 
appears that in the wide range of parameters , the zigzag AFM ordering is stable not 
only locally, but also globally, meaning that this is indeed the magnetic ground state of the 
system. We also expect the anisotropy of optical properties, which is related with the distri- 
bution of the Cg states at the corner atoms of the zigzag chain. The two-fold degeneracy of 
the Cg levels plays a very important role in the problem, and is responsible for the insulating 
behavior and the orbital ordering developed at the bridge sites of the zigzag chain. 

Thus, we believe that the nature of the low-temperature CE state in manganites is mag- 
netic. The charge ordering reduces mobility of the Cg electrons and tends to destroy the 
ferromagnetic coupling within the chain. Therefore, it is expected to be small and play only 
minor role in the problem. 

Applicability of this picture to the high-temperature regime is still in question. Particu- 
larly, many charge ordered manganites are characterized by the existence of two transition 
temperatures: the long-range AFM ordering occurs below the Neel temperature Tn, whereas 
the charge and orbital ordering orderings occur simultaneously below another critical tem- 
perature Tco- In accordance with neutron scattering measurements, it holds T^kTqo (see, 
e.g., Ref. W). This fact is typically considered as the strong experimental evidence that the 
charge and orbital degrees of freedom play a decisive role in the interval Tn<T<Tco and 
drive the AFM CE ordering at Tj<(M^ However, other experimental data suggest that the 
exact nature of the state realized between Tn and Tco is much more complicated, and still 
far from the complete understanding. 

1) Some magnetic experiments simply "do no see" the existence of the transition point Tn. 
For example, the behavior of spin magnetization in the external magnetic field shows the 
characteristic metamagnetic transition both below and above Tn, which disappears only at 
TcQ .^^i^'^ The magnetic susceptibility shows a pronounced peak at Tco and no anomaly at 
Tn.— These data suggest that, although there is no long-range magnetic order, the magnetic 
fiuctuations and the effects of short-range order continue to play a very important role in 
the interval Tn<T<Tco- 
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2) The high-resolution electron microscopy clearly reveals the existence of the one- 
dimensional chains (stripes) even above Tn, though the exact from of the stripes is not 
quite clear yet.— It is also not clear whether these stripes have a magnetic origin or not. If 
they do, it would be a strong argument in the favor of the magnetic scenario considered in 
the present work, although the perfect geometry of the FM zigzag chain may be destroyed 
at elevated temperatures. In this respect, an interesting theoretical result has been obtained 
by Hotta et al. They argued that the class of one-dimensional zigzag objects, which pre- 
serves the insulating behavior at x=0.5 is actually much wider and is not limited by the 
zigzag chain shown in Fig. ^ All these new zigzag configuration have higher energies, and 
therefore are not realized as the ground state at low temperatures. However, they may 
carry a considerable weight in the thermodynamics averages and therefore contribute to the 
high-temperature behavior of the charge ordered manganites. 
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APPENDIX A: LOCAL DENSITIES OF STATES FOR THE ZIGZAG CHAIN 

The local density of states at the site r is given by 

Ni^{e) = -hnvG^^^ie). (Al) 

TT 

Using the explicit expressions for the matrix elements of the Green function, Eqs. (13)-©, 
noting that — ^ lim^^o+ '^'oi[e—e'+irj\~^=5{e—e'), and taking into account that 

where the derivative is evaluated at the point s^=S-i one can find the following expressions 
for the partial densities of states: 

N\{e) = - , k-^cl ^^3^ 

vrv/[3 + 5(Ac-5)]K5-Ac)-l]' 



.Ti^ ^ 1 /3 + £(Ac-e) ^ , 

^lie) = . , . | J . \ . : + ^^5(6 - Ac), (A4) 
27t\s - Ac\ y e[e - Ac) - I 2 
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NU^) = . .^A + ^^'^(^ - ^c). (A5) 

The total density of states, N{e)=N\{e)+N}^{e)+N'^{e), is given by 

N{6) = -^ |2g-Ac| ^s{e-Ac). (A6) 

vr^[3 + 5(Ac-e)]Ke-Ac)-l] 

APPENDIX B: THE TOTAL ENERGY CHANGE DUE TO SMALL ROTATIONS 
OF THE SPIN MAGNETIC MOMENTS 

Here we argue that the invariance of the total energy with respect to uniform rotations of 
all spin magnetic moments by the same angle automatically leads to the isotropic Heisenberg 
form for the total energy change associated with small rotations of an arbitrary type - Eq. 
(|T9|l . For simplicity we consider the FM ordering for which e^=e^ does not depend on the 
site i. By expanding Et[{e^ + dsi}] in the Taylor series we obtain in the second order: 

E,[{e' + 6e,}] ^ E,[e'] - ^ f,{6^'^e, " e°) - ^ ^^(^(^^e, ■ ^^^^e,). (Bl) 

i ij 

In Eq. ()B1|) we employed the requirement of rotational invariance for the perturbed system: 
if all spins are rotated the same angle, the total energy does not change. Therefore, e° and 
dsi may enter the expression for the total energy change only through the scalar products. 
Noting that ^^^^e, and 6^^'^ ej are given by Eq. ()18p and employing again the requirement of 
rotational invariance for the unperturbed system (that is the total energy does not change 
if all spins {e°} are rotated by the same angle difi^=dcfi) we find the following condition: 

fi = Jij ■ 

j 

This proves validity of Eq. ()19|) for the FM ordering. The generalization to the AFM 
ordering is straightforward, because similar arguments can be applied both for the intra- 
and inter-spin-sublattice interactions in the AFM state. 
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FIG. 1: Spin, orbital and "charge" ordering in two neighboring planes of the CE-type antiferro- 

30 

magnetic structure: symbols A and A' stand for the bridge sites, symbols B and B' stand for the 
corner sites. an=\A^B^ I is the cubic lattice oarameter. 
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FIG. 2: Total and partial densities of states of the tight-binding model for Ac=0 (no charge 
disproportionation) . 
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FIG. 3: Effect of the charge disproportionation on the ferromagnetic double-exchange coupling 
within the chain. 
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FIG. 4: Total energies (bottom) and the regions of local stability (top) for various magnetic states 
in the three-dimensional (3D) and two-dimensional (2D) perovskite lattices for x=0.5. The data 
for the CE-type antiferromagnetic ordering are divided in two parts corresponding to the hard (h) 
and soft (s) regimes of the local stability. 
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FIG. 5: Nearest-neighbor interactions in the double-exchange hmit as the function of the Jahn- 
Teller distortion at the bridge sites. 




FIG. 6: Partial contributions to the conductivity tensor of the tight-binding model. 
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FIG. 7: Partial densities of Mn(3deg) (top) and Mn(4p) (bottom) states for Yi/2Sr3/2Mn04 in the 
zigzag antiferromagnetic state. Two panels show results of calculations with (right) and without 
(left) the Jahn- Teller distortion. The Fermi level (the top of occupied band) is at zero. The arrows 
t and I stand for the majority- and minority-spin states. 
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FIG. 8: Elements of the conductivity tensor for Yi/2Sr3/2Mn04 in the zigzag antiferromagnetic 
state. Different lines correspond to different values of the Jahn- Teller distortion 6. 
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